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>-very  effective  In  the  extraction  of  granite  regions  when  (l)  data  were 
in  ratio  format,  (2)  feature  variables  included  both  tone  and  texture 
information,  and  (3)  the  classifier  is  capable  of  handling  non-normal !y 
distributed  data.  Classification  errors  occurred  when  there  exists  pixel* 
of  non-g ran i te'catego ry  whose  spectral  and  texturai  properties  are 
statistically  similar  to  that  of  granite  pixels.  Two  cases  of  errors  can 
be  noted:  Type  1  pixels  located  at  the  periphery  of  the  granite  regions 
and  Type  2  pixels  located  far  away  from  the  core  of  the  granite  areas. 

To  reduce  the  error  rate,  an  unsupervised  classification  method 
based  on  the  concept  of  region  growing  and  texture  clustering  analysis 
was  employed  to  segment  the  scene  in  multiple  stages  and  thus  depict 
edge  patterns  by  the  scene  content  and  a  gradual  mathematical  generali¬ 
zation  process.  Identification  of  the  grantie  regions  becomes  a  labeling 
.  •process  using  the  training  sets  Information.  Since  the  Regions  algorithm 
is  based  on  an  additional  constraint  on  spatial  contiguity,  the  above- 
mentioned  two  types  of  errors  can  be  effectively  reduced  bacause  sharp 
edges  exist  between  the  granite  and  non-granite  pixels  in  the  study  area 
>The  final  decision  regarding  the  delineation  of  the  granite  regions 
is  based  on  the  intersection  of  two  classification  maps  using  a  simple 
map  overlay  analysis.  The  result  yields  a  correct  classification  rate 
of  about  95  percent  based  on  a  visual  comparision  between  the  composite 
classification  map  and  the  ground  truth  information  given  in  the  U.S.G.S 
geological  map  of  the  study  area,  -c ~ 

To  improve  the  developed  techniques  for  lithological  analysis,  it 
is  recommended  that  additional  experiments  be  conducted  using  other 
regions  in  the  United  States  centering  around  the  following  tasks: 

(1)  developing  algorithms  for  merging  supervised  and  unsupervised 
classification  methods; 

(2)  finetuning  the  Region  algorithm  by  adding  subroutings  to  output 
digital  information  of  each  segmented  region; 

(3)  developing  a  color  prediction  model  for  rock  types  identification 
using  the  texture  and  tone  information  in  the  color  domain  with  a 
color  monitor;  and 

(1*)  developing  change  detection  methods  for  monitoring  purposes  based 
on  the  extension  of  the  above  three  methods. 
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Executive  Summary 

It  has  been  determined  in  the  literature  on  seismology  and  geophysics 
that  the  recorded  seismic  ware  energy  from  nuclear  explosions  is  highly 
dependent  upon  the  actual  yield  of  the  explosion  and  its  interaction  with 
the  environments  in  which  the  detonation  occurs.  These  environmental  factors 
can  be  characterized  by  the  depth  of  explosion  below  the  surface,  the  degree 
of  coupling  between  the  charge  and  the  adjacent  medium,  and  the  lithological 
nature  of  the  test  sites.  Therefore,  the  analysis  of  rock  type  at  the  test 
sites  is  the  first  step  in  nuclear  monitoring. 

The  LANOSAT  data  have  been  determined  effective  for  terrain  analysis. 

The  choice  of  the  LANOSAT  imagery  for  rock  types  analysis  at  the  nuclear 
test  sites  is  also  based  upon  the  fact  that  it  can  provide  world-wide 
coverage  with  repetitive  observations  for  monitoring  purposes.  The  goal  of 
this  study  is  to  test  the  feasibility  of  utilizing  LANDSAT's  digital,  multi- 
spectral  reformation  for  rock  types  discrimination  at  the  nuclear  test  sites, 
based  on  the  texture-tone  analysis  algorithms  of  the  image  processing  systems 
at  Susquehanna  Resources  and  Environment,  Inc. 

The  experiments  were  based  on  two  subframes  of  LANDSAT  MSS  data,  cover¬ 
ing  two  geological  quadrangles  of  Nevada.  Whereas  Site  1  (Antler  Peak  Quad) 
was  used  mainly  for  methodological  development,  Site  2  (Duffer  Peak  Quad) 
was  designed  as  an  analog  test  site  to  foreign  areas  for  performing  the 
task  of  extracting  the  granite  regions. 

The  task  was  accomplished  by  using  two  separate  but  complimentary  image 
processing  techniques.  The  first  technique,  a  supervised  classification, 
was  designed  to  extract  granite  regions  using  four  ratio  bands  (4/7,  4/6, 

5/7,  and  6/7)  based  upon  four  manually  selected,  but  automatically  pre- 
processed  training  sets.  The  non-granite  regions  were  extracted  as  well 
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using  the  reject  category  of  the  classification  model.  The  second  method, 
an  unsupervised  classification  procedure  based  on  the  concepts  of  region 
growing  and  texture  analysis,  was  designed  to  delineate  granite  regions 
using  one  ratio  band  (A/7  is  most  effective)  by  growing  the  granite  regions 
from  the  cores  of  the  training  sets  to  the  edges  bordering  non-granite 
areas . 

The  final  granite  regions  were  defined  by  the  intersection  of  two 
granite  images  produced  by  two  different  image  analysis  techniques.  The 
result  indicates  that  a  very  high  level  of  correct  classification  rate-- 
95  percent  or  bettei — has  been  achieved,  based  on  an  overlay  analysis  using 
the  classification  result  against  the  geologic  map  produced  by  the  U.S. 
Geological  Survey. 

Though  the  defined  task  of  extracting  granite  regions  has  been  success¬ 
fully  accomplished,  it  is  necessary  to  test  the  developed  image  processing 
and  analysis  techniques  using  additional  U.S.  test  sites  before  they  are 
applied  to  foreign  regions.  The  reasons  are  (1)  fine  tuning  of  the  methods 
are  required  to  handle  diverse  patterns  of  lithological  associations,  and 
(2)  the  LANDSAT  imagery  can  be  exploited  further  for  detecting  environmental 
and  man-made  changes  before  and  after  nuclear  explosions. 
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Applications  of  Texture  Analysis  for  Rock  Types  Discrimination 

Section  A:  Introduction 

Ever  since  the  Soviet  Union's  detonation  of  its  first  nuclear  device 
prototypes,  both  the  realities  of  an  arms  race  and  the  requirement  to  main¬ 
tain  scientific/technological  advantages  have  forced  the  United  States  to 
expend  significant  resources  in  monitoring  of  foreign  nuclear  tests.  Sophis¬ 
ticated  technologies  that  have  evolved  about  the  framework  of  seismology  and 
geophysics  have  made  significant  contributions  in  satisfying  the  national 
requirement  to  detect,  locate,  identify  and  yi eld-quant i fy  world-wide  nuclear 
detonations.  Yet,  there  is  room  for  improvement  using  non-seismic  methods, 
particularly  in  the  area  of  yield  estimation.  To  this  end,  this  study  is 
intended  to  develop  image  processing  and  analysis  methodologies  for  the 
discrimination  and  identification  of  rock  types  at  nuclear  test  sites.  The 
rationaleof  this  approach  is  based  on  the  fact  that  the  recorded  seismic 
wave  energy  resulted  from  nuclear  explosion  depends  on  the  following  environ¬ 
mental/lithological  factors: 

(1)  the  actual  yield  of  the  explosion; 

(2)  depth  of  the  explosion  below  the  surface; 

(3)  the  degree  of  physical  coupling  between  the  charge  and  the 
adjacent  medium;  and 

(*0  the  geological  nature  of  the  median  in  which  the  detonation 
occurs . 

Indeed,  rock  types  analysis  is  the  first  step  in  yield  estimation. 

To  accomplish  the  goal  of  rock  types  discriminated  at  the  nuclear  test 
sites,  LANDSAT's  mul t i spectra  1  data  were  used.  The  choice  of  the  LANDSAT 
imagery  is  based  on  the  fact  that  it  is  capable  of  providing  a  world-wide 
and  repetitive  coverages  and  thus  a  basis  for  monitoring  nuclear  test 
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activities.  The  thrust  of  this  study  is  to  exploit  the  digital  information 
of  LANDSAT  data  in  the  context  of  texture-tone  analysis  for  such  purposes. 

To  test  the  feasibility  of  the  SR&E's  image  processing  system  for  litho¬ 
logical  analysis,  two  test  sites  in  Nevada  were  utilized,  the  Antler  Peak 
Quadrangle,  Nevada  at  the  scale  of  1:62,500  (ANA1)  and  the  Duffer  Peak 
Quadrangle,  Nevada  at  1:48,000  (ANA2)  as  analogs  to  foreign  nuclear  test 
sites.  Specifically,  the  first  site  (ANA1)  was  used  as  a  testbed  for  methodo¬ 
logical  development;  whereas  the  second  site  (ANA2)  was  designed  as  an  analog 
area  for  extracting  granite  regions. 

To  classify  granite  versus  non-granite  regions,  two  complementary  image 
analysis  techniques  were  employed.  First,  a  supervised  classification  analy¬ 
sis  was  conducted  to  delineate  granite  areas  based  on  manually  selected,  but 
digitally  pre-processed  training  sets,  and  to  reject  non-granite  regions 
based  on  a  pre-set  statistical  model/probability  level  for  identifying  pixels 
which  are  significantly  different  from  the  training  sets.  Second,  an  un¬ 
supervised  clustering  analysis  based  on  SRSE's  Region  Growing  Texture  Clus¬ 
tering  algorithm  was  performed  to  extract  granite  areas  by  region  growing 
from  the  cores  of  the  granite  training  sets.  The  final  definition  of  granite 
regions  is  based  on  the  intersect  ion  of  these  two  sets  of  "gran  i  te  maps." 

Section  B:  A  Brief  Review  of  Relevant  Literature 

Prior  to  1972  and  the  launch  of  LANDSAT,  pioneer  work  on  reflective 
properties  of  minerals  was  accomplished  by  Hunt  and  Salisbury  at  the  USAF 
Cambridge  Research  Laboratories  (1970,  1973)-  Their  study  and  explanation 
of  reflective/transmission  properties  of  both  minerals  and  rocks  in  the  visible 
and  near- inf rared  regions  serves  as  a  basis  for  semi-automatic  rock  discrimi¬ 
nation  to  techniques  that  exploit  the  spectral  (tone)  parameters  of  multi- 
spectral  imagery.  Rather  than  being  a  simple  empirical  result,  it  turns  out 


that  minerals  and  rocks  spectral  characteristics  are  a  direct  function  of 
the  physics  and  theory  associated  with  crystal-field  theory  (Burns,  1970), 
as  evidenced  from  theoretical  and  laboratory  analyses  of  the  rocks  and  minerals 
of  the  moon  (McCord,  1968;  further  McCord,  et  a  1 ,  1972). 

Since  then,  scientists  at  the  U.S.  Geological  Survey,  Jet  Propulsion 
Laboratory  and  NASA/Goddard  Space  Flight  Center,  have  attempted  to  exploit 
LANDSAT  MSS  data  for  rock  types  analysis  as  evidenced  from  Goetz,  et  al 
(1973),  Goetz,  et  al  (1975),  Vincent, et  al  (1975),  Rown,  et  al  (1976), 

Rown,  et  al  (1977)/  Abrams,  et  al  (1977),  and  Podwysocki ,  et  al  (1977). 

Recent  work  by  other  researchers  including  Lyon  (1977),  Lyon,  et  al  (1978), 
Hunt  (1977),  and  Siegrist  et  al  (1980),  also  emphasized  digital  processing 
of  LANDSAT  and  other  types  of  mul tispectral  scanner  data  for  optimal  combi¬ 
nation  of  spectral  channels  for  rock  discrimination. 

While  the  majority  of  the  work  cited  above  emphasized  rock  types  analy¬ 
sis  and  identification  with  color  enhancement  techniques  with  LANDSAT  images, 
our  study  is  devoted  exclusively  to  extracting  rock  types  using  the  digital 
information  of  the  LANDSAT  MSS  data  in  the  context  of  texture  analysis,  which 
has  been  largely  neglected  by  previous  researchers. 

Section  C:  The  Tasks,  Data  Analysis  and  Results 
1 .  Tasks  to  be  Accomplished 

a.  Preliminary  Testing  on  the  Proposed  Methodologies 

At  the  beginning  of  our  research,  the  task  was  loosely  defined 
as  discrimination  of  rock  types  with  LANDSAT  data  using  our  texture 
analysis  algorithms. 

Using  more  than  20  training  sets  evenly  distributed  over  the 


entire  Site  1  (Antler  Peak  Quad.),  it  was  determined  that  (1)  our 
texture  algorithm  is  capable  of  separating  these  training  sets  with 
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a  correct  classification  rate  of  over  95  percent  even  though  many 
of  the  training  sets  belong  to  the  same  major  rock  types--sed imentary , 
igneous  and  metamorphic;  and  (2)  while  the  three  major  rock  types  are 
well  separated,  in  terms  of  statistical  means  of  texture-tone  variables, 
the  range  of  these  measurements  from  training  sets  within  the  same 
major  rock-type  is  not  small  for  all  training  sets,  meaning  that  local 
variations  exist  in  major  rock  types. 

In  the  mapping  analysis,  it  was  determined  that  our  region-growing 
texture  clustering  algorithm  is  apparently  effective  in  delineating 
surface  material  which  may  be  relatable  to  the  bedrock  information  with 
LANDSAT  data  without  compression.  It  is  not  effective  when  the  LANDSAT 
data  are  compressed  by  a  factor  of  three, 
b.  Tasks  Determined  for  the  Phase-1  Effort 

From  a  discussion  session  among  Mr.  Best  of  AFOSR,  Col.  Lowrey  of 
DARPA,  Mr .  Rachl in  and  his  colleagues  of  U.S.  Geological  Survey,  and 
Or.  Hsu,  it  was  determined  that  the  tasks  of  our  effort  should  be  aimed 
at  answering  the  following  three  questions: 

1.  How  well  can  we  map  the  granite  areas  versus  non-granite 
regions  using  our  supervised  and  unsupervised  classification 
methods  in  the  context  of  texture  analysis? 

2.  What  are  the  factors  affecting  the  classification  results-- 
slope,  drainage  pattern,  data  used,  methodologies  utilized? 

3.  What  are  the  potential  contribution  of  image  processing 
techniques  and  methodologies  towards  the  discrimination  and 
even  identification  of  rock  types  using  LANDSAT  data? 

For  data  analysis,  a  study  area  within  Duffer  Peak  Quadrangle,  Nevada 
was  selected  by  Mr.  Dempsie  of  U.S.  Geological  Survey.  Furthermore,  based 
on  the  geologic  map,  22  training  sets  were  selected  manually  to  cover  four 


major  rock  types:  (1)  granite,  (2)  metamorphic,  (3)  volcanic,  and  (4) 
unconsol i dated . 
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2.  The  Data  Sets 


Corresponding  to  the  study  area  selected  by  Mr.  Dempsie  of  U.S.G.S., 
a  digital  set  composed  of  (256  x  256)  pixels  of  LANDSAT  MSS  data  was 
determined  using  visual  analysis.  The  Northwest  corner  of  the  data 
set  is  located  at  (row  1268,  column  1987). 

In  addition,  the  training  sets  with  their  locations  within  the 
(256  x  256)  frame  have  also  been  determined  as  follows: 

1 .  Group  1 :  Gran i te 

Gl:  (15,124),  (15,144),  (27,144),  (27,124) 

G2:  (56,  38),  (56,  55) ,  (71 ,  55) ,  (71  ,38) 

G3:  (130,133),  030,147),  (142,147),  (142,133) 

G4 :  (132,7),  (132,19),  (146,19),  (146,7) 

G5:  (170,137),  (170,150),  (181,150),  (181,137) 

2.  Group  2:  Metamorphic 

Bl:  (75,158),  (75,172),  (85,172),  (85,158) 

B2 :  (115,163),  (115,182),  (129,182),  (129,163) 

B3:  (150,175),  (150,189),  (166,189),  (166,175) 

B4:  (95,89),  (95,106),  (105,106),  (105,89) 

3.  Group  3:  Volcanic 

HI:  (213,144),  (213,159),  (221,159),  (221,144) 

H2 :  (215,174),  (215,190),  (223,190),  (223,174) 

H3:  (241,186),  (241,197),  (248,197),  (248,186) 

H4:  (217,13),  (217,25),  (228,25),  (228,13) 


FI:  (188,81),  (188,93);  (197,93),  (197,81) 


F2:  (228,106),  (228,117),  (235,117),  (235,106) 
F3:  (213,124),  (213,134),  (222,134),  (222,124) 
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Tl:  (234,68),  (234,80),  (243,80),  (243,68) 

T2:  (246,100),  (246,112),  (254,112),  (254,100) 

4.  Group  4:  Unconsolidated 

Qol:  (132,51),  (132,64),  (139,64),  (139,51) 

Qfl:  (51,196),  (51,207),  (61,207),  (61,196) 

Qf2 :  (150,212),  (150,223),  (158,223),  (158,212) 

Q.1 :  (56,237),  (56,249),  (65,249),  (65,237) 

b.  Derived  Data  Sets  to  be  Analyzed 

To  remove  the  shadow  effect  of  the  original  LANDSAT  data,  and  to 
extract  information  from  four  MSS  bands  simultaneously,  the  following 
date  sets  are  generated. 

(1)  First,  second  and  third  components  from  the  four  MSS  bands; 

(2)  Six  ratio  bands  from  the  four  MSS  bands:  4/5,  5/6,  6/7, 

4/6,  4/7,  and  5/7. 

(3)  The  first  component  map  from  4  selected  ratio  bands. 

Therefore,  ten  derived  image  data  sets  are  available  for  analysis  in 

addition  to  the  original  four  MSS  bands.  The  location  of  the  training  sets 
with  respect  to  these  derived  data  sets  remain  the  same. 

3 •  Image  Processing  and  Data  Analysis  Methodologies  Utilized 

To  analyze  the  relationship  between  the  selected  training  sets,  and  to 
classify  the  granite  areas  versus  non-granite  regions,  the  following  analyti¬ 
cal  techniques  are  utilized. 

a .  Extraction  of  texture-tone  information  of  the  training  sets  and 
the  enti re  data  set. 

Using  the  texture-tone  extraction  algorithm,  23  texture-tone- 
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ratio  variables  have  been  generated  for  any  given  pixel  from  four 
mul ti -spectral  bands  using  (3  x  3)  moving  grid.  They  are  composed 
of  k  tone  variables,  12  texture  variables  (3  from  each  band),  1 
linear  feature  variable,  and  final  6  ratio  variables. 

For  data  analysis,  the  analyst  is  able  to  select  a  portion  of 
the  23  variables. 

b.  Analysis  of  the  Training  Sets 
Based  upon  the  selected  variables  from  23-variable  system,  typi¬ 
cally  we  use  three  variables,  the  training  sets  will  be  analyzed  and 
edited  so  that  each  training  set  will  meet  the  following  two  criteria: 

(1)  single  mode;  if  two  modes  exist  in  one  training  set,  the  set 
will  be  split  into  subsets; 

(2)  extreme  outliers  are  to  be  removed  based  on  a  statistical 
confidence  level. 

c.  Discriminant  Analysis  of  the  Training  Sets 
After  the  training  sets  are  edited  or  preprocessed,  they  will  be 

analyzed  in  terms  of  how  close  they  are  between  pairs  of  training  sets 

using  the  means  of  selected  tone-texture  variables.  The  distance  is 

2 

generally  measured  by  statistical  distance  called  Mahalanobis  D  with 

or  without  a  log-determinant  term. 

2 

While  the  D  distance  is  indicative  of  the  degree  of  similarity  and 
dissimilarity  between  two  training  sets,  the  analyst  usually  uses  a  con¬ 
fusion  matrix — classification  result  using  only  the  training  sets--to 
examine  how  well  these  training  sets  are  separated.  The  analyst  will 
then  decide  whether  he  should  proceed  with  a  classification  analysis 
of  the  entire  test  set.  In  general,  if  dissimilar  training  sets  are 


confused,  a  classification  analysis  should  not  be  conducted. 
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d.  Supervised  Class! ficatlon  Methods 

As  mentioned  earlier,  a  supervised  classification  analysis  can 
be  made  only  when  the  training  sets  are  well  separated.  To  achieve 
this  goal,  the  following  steps  can  be  taken: 

(1)  purify  the  training  sets  as  in  (b) ; 

(2)  change  the  location  of  the  training  sets; 

(3)  increase  the  power  of  the  feature  extractor  by  using 
(i)  more  texture-tone  variables,  and  (ii)  using  dif¬ 
ferent  spectral -band  combinations;  and 

(4)  increase  the  power  of  the  classifier  by  using  a  non- 
Gauseian  model  if  the  data  are  essentially  non-multi¬ 
variate  normal . 

In  the  analysis,  we  have  done  all  thse  image  processing  techniques  except 
step  (2),  changing  the  location  of  the  training  sets. 

e.  Scene  Segmentation  with  Unsupervised  Training/Classification  Method 
This  analysis  is  intended  to  extract  the  granite  regions  first 

based  on  segmentation  concept  using  a  region-growing  texture  clustering 
algorithm.  Once  the  entire  test  area  is  segmented  into  numerous  sub- 
regions  according  to  different  1 evel s  of  threshold i ng  (generalization), 
we  are  able  to  extract  the  granite  regions  according  to  the  location 
of  the  training  sets. 

Since  this  algorithm  is  based  on  local  statistics  or  edge  informa¬ 
tion  instead  of  global  separation  employed  by  the  supervised  classifi¬ 
cation  method,  it  should  be  used  as  a  complementary  classification 
method  instead  of  a  replacement  of  the  supervised  method. 

f .  Comparative  Analysis  with  the  Results  from  Supervised  and  Unsuper- 
vised  Classification  Method 
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g.  Analysis  of  Factors  Influencing  the  Classification  Results 

In  this  analysis  we  will  concentrate  our  effort  on  two  broad 
categories:  (1)  classification  results  influenced  by  terrain  factors, 

surfacial  material,  slope,  drainage,  etc.,  and  (2)  classification  re¬ 
sults  influenced  by  the  data  sets  and  techniques,  feature  extractors, 
and  classifiers,  used  in  the  analysis. 

4.  Experiments  Conducted  and  Research  Results 

The  description  of  this  section  regarding  data  analysis  corresponds  to 
the  methodologies  discussed  in  the  previous  section. 

a.  Generation  of  Texture-tone  Variables  for  Data  Analysis  by  Super¬ 
vised  Classification  Methods 

Using  the  original  LANDSAT  MSS  data  and  dervied  ratio  bands, 
these  data  sets  containing  23  texture-tone  variables  for  each  pixels 
were  generated: 

(1)  Data  Set  1  is  composed  of  MSS  Bands  4,  5,  6,  and  7; 

(2)  Data  Set  2  is  composed  of  4  ratio  bands:  4/5,  5/6,  6/7  and  4/7;  and 

(3)  Data  Set  3  is  composed  of  4  ratio  bands:  4/7,  4/6,  5/7  and  6/7. 

b.  Analyses  of  the  Training  Sets 

The  training  sets  selected  manually  by  Mr.  Dempsie  of  U.S.  Geologi¬ 
cal  Survey  were  analyzed  for  two  distinctive  purposes: 

(1)  All  22  training  sets  were  analyzed  to  detect  the  confusion 
pattern  between  the  granite  sets  and  the  non-granite  sets;  and 

(2)  Four  granite  sets  of  the  original  five  granite  sets  were  pre- 
processed  for  serving  as  calibration  samples  for  classifica¬ 
tion  analysis. 


i 


The  aniaysis  starts  with  preprocessing  of  the  training  sets  aimed  at 
detecting  whether  bi-modal  distribution  and  outliers  exists  in  each  set. 
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Indeed,  it  was  determined  that  in  the  original  MSS  data  sets  B2,  B3,  B4, 

T1  and  QqI  are  bi-modal,  and  thus  each  was  split  into  two  subsets,  resulting 
in  28  training  sets,  instead  of  22  sets  in  original  design.  Furthermore, 
outliers  are  edited  and  eliminated  using  a  predetermined  distribution  scheme. 
Table  1  summarizes  the  training  sets  information,  original  and  after  pre¬ 
processing,  whereas  Appendix  1  gives  the  sample  statistics  of  the  texture- 
tone  variables  for  each  training  set. 

Using  the  statistics  given  in  Appendix  1,  an  analysis  was  conducted  to 
reveal  the  confusion  pattern  between  the  granite  training  sets  (set  1  through 
set  5)  and  the  rest.  The  results  are  given  in  Table  2,  and  it  indicates  that 
(1)  granite  set  1  is  highly  confused  with  non-granite  sets,  and  (2)  granite 
set  2  through  set  5  are  highly  correlated  among  themselves,  but  are  not  con¬ 
fused  with  non-granite  sets  (set  6  through  set  28).  For  instance,  the  per¬ 
centage  of  correct  classification  of  G1  into  G1  through  G5  is  only  61.7, 
whereas  the  figures  for  G2,  G3,  G4  and  G5  are  100,  86.1,  99*5  and  99-4, 
respectively.  It  was  therefore  determined  that  training  set  G1  should  be 
eliminated  from  the  design  sets  in  the  final  classification  analysis  aimed 
at  delineating  granite  versus  non-granite  regions. 

The  same  analyses  were  also  applied  to  the  training  sets  with  ratio 
data,  instead  of  the  original  MSS  data.  Whereas  Appendix  2  gives  the  texture- 
tone  statistics  of  each  training  set.  Table  3  summarizes  the  result  of  pre¬ 
processing;  from  the  original  22  sets,  25  sets  are  obtained,  instead  of  28 
as  in  the  case  of  the  original  MSS  data.  This  means  that  the  ratio  data 
are  more  homogeneous  than  the  raw  MSS  data  because  the  LANDSAT's  shadow 
effect  has  been  removed  by  the  ratioing  process. 

In  terms  of  confusion  pattern  between  granite  sets  versus  non-granite 
sets,  no  significant  difference  exists  between  the  raw  MSS  data  and  the 
ratio  data,  as  indicated  in  Table  4;  namely,  granite  set  1  is  totally  confused 
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with  other  non-granite  sets,  granite  set  2  through  granite  set  5  are  similar 
among  themselves,  but  quite  different  from  other  non-granite  rocks. 


Raw  Data 

Table  1:  Preprocessing  of  the  Training  Sets  with  Original  MSS  Data 

(Band  4  and  Band  7) 


Original 

Training  at 
ID  # 


Training 
*'0c^e  it  of  Points 


Set  Size 

(2) 

after 

Preprocess i ng 


New  Set  ID 

New  Traini ng 
Set 
ID  a 


1 

Granite  G1 

273 

269 

1 

2 

G2 

288 

z78 

2 

3 

G3 

195 

190 

3 

4 

G4 

195 

187 

4 

5 

G5 

168 

164 

5 

6  Metamorphic  B1 

165 

162 

6 

7 

B2 

o 

o 

<*\ 

141 

154 

7 

8 

8 

B3 

165 

95 

70 

9 

10 

9 

B4 

198 

91 

92 

11 

12 

10 

Volcanic  HI 

144 

143 

13 

11 

H2 

153 

150 

14 

12 

H3 

96 

96 

15 

13 

H4 

156 

93 

63 

16 

17 

14 

FI 

130 

126 

18 

15 

F2 

96 

94 

19 

16 

F3 

no 

108 

20 

T1 

130 

64 

21 

17 

51 

22 

18 

T2 

117 

111 

23 

Unconsol idated 

19 

Qfi 

132 

131 

24 

20 

Qf2 

108 

105 

25 

21 

Q1 

130 

130 

26 

22 

V 

112 

58 

50 

27 

28 

i 
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Table  2:  Confusion  Matrix:  Percent  of  Correct  Classification 

With  the  Raw  Data 


New 

10  # 

1 

2 

3 

4 

5 

Gi 

6.06 

0 

1.1 

0 

0 

G2 

0 

82.4 

0 

17.3 

0.4 

G3 

0 

0 

84.7 

0.5 

9-0 

G4 

0 

18.2 

0.5 

75.4 

3.7 

G5 

0 

4.9 

6.7 

0 

81.1 

10 

11 

12 

13 

14 

2.6 

0 

0.4 

l.l 

8.2 

0 

0 

0 

0 

0 

0 

0 

1.6 

0 

0.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.6 

New 

ID  # 

15 

16 

17 

18 

19 

20 

Gl 

1.9 

3.4 

0 

0.4 

1.5 

0.4 

G2 

0 

0 

0 

0 

0 

0 

G3 

1.1 

0.5 

0 

1.1 

0 

0.5 

G4 

0 

0 

0 

0 

0 

0 

G5 

0 

0 

0 

0 

0 

0 

21  22  23  24  25  26  27  28 

0  0.4  2.2  0  2.2  2.2  0  0 

00  000  0  00 

0  0  0  0  0  0  0  0.5 

0  0  0  0  0  0  0  0.5 

00  000  0  00 


Table  3:  Preprocessing  of  Training  Sets  with  Ratio  Data 


Training  Set  10 
If 

1  Granite 

Code 

H  of 
Points 

After  pre¬ 
processing 

New 
ID  # 

G1 

273 

258 

1 

2 

G2 

288 

278 

2 

3 

G3 

195 

186 

3 

4 

G4 

195 

178 

4 

5 

G5 

168 

162 

5 

6  Metamorphic 

8! 

165 

165 

6 

7 

B2 

300 

I83 

111 

7 

8 

8 

B3 

165 

74 

83 

9 

10 

9 

84 

198 

196 

11 

10  Volcanic 

HI 

144 

143 

12 

11 

H2 

153 

150 

13 

12 

H3 

96 

93 

14 

13 

H4 

156 

154 

15 

14 

FI 

130 

127 

16 

15 

F2 

96 

93 

17 

16 

F3 

110 

108 

18 

17 

T1 

130 

127 

19 

18 

T2 

117 

115 

20 

19  Unconsol  idated  Qf  I 

132 

130 

21 

20 

Qf2 

108 

105 

22 

21 

Q1 

130 

63 

57 

23 

24 

22 

Qo> 

112 

1 1 1 

25 

14 


Table  4:  Confusion  Matrix  Between  Granite  and  Non-Granite  Sets  with 
Ratio  Data  (B4/B7,  B4/B6,  B5/B7,  B6/B7) 


1  2  3  4  5-  6  -7  8  9  10  11  12 


G1 

13.2 

0 

0 

0 

1.2 

0.4 

2.7 

1.9  1.6 

0.4 

13.2 

1.2 

G2 

0 

48.2 

0 

13.7 

25.5 

0 

0 

0 

0 

0 

0.4 

0 

G3 

0 

0 

32.8 

0 

61.3 

0 

0 

0 

0 

0 

1.6 

0 

G4 

0 

7.3 

0.6 

37.1 

31.5 

0 

0 

0 

0 

0 

0 

0 

G5 

0 

0 

0.6 

0 

99.4 

0 

0 

0 

0 

0 

0 

0 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

G1 

3.5 

3.5 

0.4 

6.2 

7.0 

5.4 

8.1 

27-9 

0 

1 .6 

0.4 

0.4 

0 

G2 

0 

0 

0 

0 

0 

0 

0 

0 

2.5 

0 

0 

0 

9.7 

G3 

9.5 

2.7 

0.5 

0 

0 

0.5 

0 

0 

0 

D 

0 

0 

0 

G4 

0 

0 

0 

0 

0 

0 

0 

0 

0.6 

EL 

0 

0 

23-0 

G5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

In  terms  of  the  confusion  pattern  among  these  4  granite  sets,  it  was 
determined  that  they  are  related  only  to  a  certain  degree  since  tbe  correct 
classification  rate  with  the  form  sets  reach  at  a  high  level  of  75  percent 
(Table  5).  This  means  significant  local  variation  of  granite  exists  in  the 
study  area. 

Table  5-  Confusion  Matrix  of  the  Granite  Training  Sets 


G2 

G3 

G4 

G4 

G2 

211 

0 

63 

5 

G3 

0 

164 

2 

20 

G4 

62 

2 

113 

2 

G5 

12 

31 

19 

102 

Correct 

Classi f ier 

CO 

r*- 

II 

.02% 

c.  Classification  Analyses  with  a  Supervised  Training  Approach 

The  goals  of  this  anlayses  are  first  to  generate  classification 
maps  of  granite  versus  non-granite  regions  using  different  data  sets, 
different  feature  variables  and  different  classifiers,  and  second  to 
compare  these  classification  results  against  the  ground  truth  infor¬ 
mation  including  terrain  information  and  bedrock  geologic  map. 

For  the  study,  numerous  experiments  have  been  conducted;  the 


following  table  indicates  the  experiments  derived  from  combinations 
of  different  data  sets  with  different  methodologies. 


Table  6. 


Experiments  of  Supervised  Classification 


Data  Sets 

Data  Sets  Data  Set  1 

Data  Set  2 

Data  Set  3 

(4  original  MSS) 

Methodology 

(4  ratio  bands: 
4/5, 4/5, 5/7, 4/7) 

(4  ratio  bands: 
4/7, 4/6, 5/7. 6/7) 

1.  Gaussian  Classi-  Exp  1  with  5 
f ier  wi th 7  tex-  granite  sets 

ture-tone  Exp  2  with  5 

variables  granite  sets 

plus  5  auto¬ 
matically  se¬ 
lected  new  sets 

Exp  1  with  4 
granite  sets 

Exp  2  with  4 
granits  sets 
plus  2  auto¬ 
matical  ly  se¬ 
lected  new  sets 

Exp  1  with  4 
granite  sets 

2.  Gaussian  Classi¬ 
fier  with  16  tex¬ 
ture  variables 

Exp  1  with  4 
granite  sets 

3.  Non-Gaussian  clas- 
s i f iers  wi th  7 
texture- tone  variables 

Exp  1  with  4 
granite  sets 

Exp  1  with  4 
granites  sets 

From  these  experiments,  it  can  be  concluded  that: 

(1)  For  rock  type  analysis.  Data  Set  3  composed  of  these  4  ratio  bands: 
4/6,  4/7,  5/7,  6/7  is  most  effective.  Data  Set  1  with  the  4  original 
4  LANDSAT  MSS  bands  is  least  affective.  Figure  1,  a  decision  map, 
indicates  that  the  vast  majority  of  granite  areas  are  correctly 
identified,  except 

(a)  the  granite  area,  within  which  the  training  set  Gl 
(which  was  not  used  in  the  analysis  is  located,  is 
largely  classified  as  non-granite,  and 

(b)  one  "metamorphic  rock"  area  as  labeled  in  the  geologic 
map  was  largely  classified  as  "granite." 

These  two  regions  will  be  investigated  further  using  our 
unsupervised  segmentation  algorithm  in  the  next  section. 

(2)  Regarding  the  classifiers,  our  Non-Gauusian  classifier  with  7 
texture-tone  variables  is  superior  to  the  Gaussian  classifier 
no  matter  whether  it  utilizes  7  or  16  texture-tone  variables. 


(3)  There  is  little  difference  between  7_variable  Guassian  classi¬ 
fier  and  16-variable  Gaussian  classifier  in  terms  of  the  con¬ 


fusion  matrix  using  the  training  sets  data. 

(4)  In  terms  of  correct  classification  of  the  granite  versus  non¬ 
granite  regions  (areal  distribution),  both  our  Gaussian  and  non- 
Gaussian  classifiers  achieved  a  level  of  over  90  percent  hit-rate. 
The  Non-Gaussian  Classifier  is  slightly  better  than  the  Gaussian 
Classifier  in  these  experiments, 
d.  Feature  Extraction  with  an  Unsupervised  Training  Approach 


The  goal  of  these  analyses  is  to  extract  homogeneous  regions  in 
the  study  area  from  various  LANDSAT  ratio  bands  using  our  region¬ 
growing  texture  clustering  analysis  algorithm.  Identif ication  of 
the  granite  regions  becomes  a  labeling  process  using  training  sets 
information  and  other  related  statistical  and  terrain  characteristics 
data.  It  is  our  intention  to  use  the  results  from  this  unsupervised 
classification  method  to  investigate  the  areas  of  misclassif ication 
by  the  supervised  classifier. 

From  the  four  LANDSAT  MSS  bands,  it  is  possible  to  derive  six 
ratio  bands:  4/5,  4/6,  4/7,  5/(?,  4/7,  and  6/7.  In  our  earlier 
experiments,  it  was  determined  that  these  ratio  bands  are  not  effec¬ 
tive  for  segmenting  regions  of  the  study  area:  4/5,  5/6  and  6/7. 

Hence,  our  experiments  for  rock-type/surface  material  analyses  utilized 
information  from  these  three  ratio  bands:  4/7,  5/7,  and  4/6.  These 
three  ratio  bands  in  fact  characterize  the  contrast  between  the  visible 
and  the  infrared  spectrum. 

To  extract  rock-type/surface  material  regions,  the  following 
experiments  were  conducted  using  our  unsupervised  classification 
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("Region")  algorithm: 


Experiment 

Data  Set 

Thresholding  Parameter 

No.  of  Passes 

1 

A2R47F3 

(2,2) 

through  (2,5) 

5 

2 

A2R47F3 

(3,3) 

through  (3,6) 

4 

3 

A2RA7F3L 

(2,2) 

through  (2.6) 

5 

4 

A2R57F3 

(3,3) 

through  (3,7) 

5 

5 

A2R57F3L 

(2,2) 

through  (2,5) 

4 

6 

A2R46F3 

(2,2) 

through  (2,7) 

7 

7 

A2R46F3L 

(2,2) 

through  (2,4) 

3 

For  the  data  set  ID, 

A2  means  Analog  Area  2 

of  the  study  areas 

R47  means 

ratio  between 

Band  4 

and  Band  7;  and  F3  means  the  third 

data  file 

(Analog  2)  and 

L  i  n 

the  data  set  name 

identifies  the 

fact  that 

the  data  set  i 

s  derived  from  double  log  ratio  mode. 

Regarding  the  thresholding  parameters,  the  first  parameter 
stands  for  the  first  stage  cutoff  regarding  the  difference  between 
adjacent  (pairs)  pixels  or  clusters;  whereas  the  second  parameter 
refers  to  the  second  stage  cutoff  for  grouping  clusters  in  terms 
of  a  weighted  geometric  distance  computed  from  a  tone  and  a  tex¬ 
ture  variable. 

Since  we  design  the  algorithm  to  perform  a  dynamic  cueing  task, 
the  analytical  results  can  be  printed  (output)  at  any  given  stage 
of  clustering  process  according  to  the  second  stage  cut-off  para¬ 
meter  specified  by  the  analyst.  Furthermore,  the  statistics  for 
each  region  of  a  given  pass  can  be  extracted  and  displayed. 

(2)  The  Results  of  the  Analyses 

By  examining  the  results  from  these  experiments,  it  can  be 
concluded  that: 

(a)  General  conclusions: 

(i)  In  general,  the  "Region"  algorithm  is  able  to  ex¬ 
tract  spatially  contiguous,  homogeneous  regions  of 
rock-type/surface  material  regions  as  defined  by 
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the  contrast  between  visible  and  infrared  spectrum 
of  the  LANOSAT  data; 

(ii)  The  performance  of  a  given  ratio  band  is  not  uniform 
over  the  entire  study  area,  meaning  that  it  may  take 
two  or  more  ratio  bands  to  extract  all  of  the  "dis¬ 
tinctive"  regions  in  the  study  area. 

(iii)  Our  "dynamic  cueing"  approach  is  able  to  reveal  the 
strength  and  weakness  of  the  contract  lines  or  zones 
between  two  adjacent  rock-type/surface-material  re¬ 
gions.  This  means  that  the  existence  of  "contact 
lines  or  zones"  is  both  spectral/spatial  information 
dependence,  and  thresholding  parameter  dependence 
as  wel 1 . 

(iv)  In  the  area  where  the  supervised  training  classifier 
failed  to  identify  the  granite  and  non-granite  re¬ 
gions,  the  unsupervised,  region-growing  algorithm 
is  capable  of  identifying  them  as  distinctive  regions, 

(v)  Combining  the  results  from  the  supervised  and  the 
unsupervised  classification  approaches,  we  believe 
that  the  correct  classification  rate  of  granite 
versus  non-granite  region  is  about  95  percent. 

(vi)  Therefore,  it  can  be  concluded  that  our  two  classi¬ 
fication  algorithms  are  indeed  complementary  for 
extracting  distinctive  rock-type/surface-material 
regions. 

(b)  Conclusions  from  individual  experiments: 

(i)  Results  from  Log  Ratio  of  Band  4  and  B7  (Figure  2). 

The  base  map  of  Figure  2  is  the  Pass  3  results  of 
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the  unsupervised  classification  with  log  ratio  of 
Band  4  and  Band  7.  Using  the  location  of  the  train¬ 
ing  sets  G1  through  G5,  the  major  granite  areas  are 
identified  and  colored  in  light  red,  whereas  the  bed¬ 
rock  regions  of  granite  are  shaded  in  green.  In 
general,  there  is  a  high  degree  of  agreement  between 
the  segmented  regions  and  the  bedrock  boundaries. 
Particularly,  by  comparing  Figure  2  against  Figure  1, 
we  are  able  to  derive  that  (1)  the  rejected  granite 
G1  area  can  be  delineated  by  the  Region  algorithm, 
and  (2)  the  confused  area  in  Figure  1  between  G2  and 
G3  can  be  discriminated  as  well.  Similar  to  the 
supervised  classifier,  the  Region  algorithm  failed  to 
distinguish  the  bedrock  granite  from  the  surfacia 1 
granite  in  the  area  near  training  set  G4.  As  will  be 
noted  later,  this  boundary  is  detected  in  the  analysis 
with  the  data  set  of  log  ratio  of  Band  4/Band  6. 

(ii)  Results  from  Log  Ratio  of  Band  4  and  Band  7  with  a 
Larger  First  Stage  Cutoff.  This  experiment  was  in¬ 
tended  to  reveal  the  effect  of  using  a  larger  first 
stage  cutoff  parameter  as  compared  to  the  above  experi¬ 
ment.  The  results  indicate  that  with  a  larger  cutoff 
the  algorithm  failed  to  detect  the  boundary  between 
granite  versus  non-granite  in  the  area  where  training 
set  G1  is  located;  the  rest  of  the  results  remains 
essentially  the  same. 

(iii)  Results  from  Double  Log  Ratio  of  Band  4  and  Band  7 
(Figure  4).  This  experiment  shows  the  effect  of 
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using  a  double  log  transformation  of  the  data  set 
Band  4/Band  7.  By  comparing  the  result  against  Figure 
2,  it  can  be  noted  that  this  double  transformation 
in  fact  has  less  discrimination  power  as  evidenced 
from  the  fact  that  there  is  boundary  shift  in  the 
area  near  training  sets  G2  and  G5. 

(iv)  Results  from  Log  Ratio  of  Band  4  and  Band  6.  In  this 
experiment,  we  replaced  Band  7  with  Band  6  of  the 
LANDSAT  infrared  spectra  and  kept  Band  4  as  a  con¬ 
stant.  The  result  shows  that  (1)  log  ratio  of  Band  4 
and  Band  6  is  able  to  detect  the  boundary  between 
bedrock  granite  and  surfacial  granite  as  evidenced 
from  the  area  near  the  location  of  G4,  but  (2)  it 
failed  to  detect  the  granite  versus  non-granite 
boundaries  in  the  areas  where  G1  and  G4  are  located. 

In  general,  it  appears  that  the  spectral  data  from 
ratio  of  Band  4/Band  6  contain  a  component  which  is 
affected  by  the  drainage  pattern  of  the  area. 

(v)  Results  from  Log  Ratio  of  Band  5  and  Band  7. 

This  experiment  shows  the  effect  regarding  a  change 
in  the  visible  band  in  the  analysis,  and  the  results 
indicate  that  the  ratio  of  Band  5  and  Band  7  is  less 
effective  for  rock  type  analysis  as  compared  to  B4/B7. 
It  appears  that  the  B5/B7  data  set  contains  a  component 
which  is  highly  affected  by  the  drainage  patterns  of 
the  area. 

(vi)  Results  from  Double  Log  Ratio  of  Band  5  and  7. 

In  general  this  experiment  indicates  that  in  certain 
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Section  0:  Image  Processing  Techniques 
Towards  Identification  of  Rock  Types 

l.  Integration  of  Supervised  and  Unsupervised  Classification  Methods. 

As  indicated  in  Section  C,  the  task  of  identifying  rock  types  can  be 
achieved  provided  that  ground  truth  information  of  the  training  areas  is 
known.  The  task  can  be  achieved  by  using  either  a  supervised  or  unsuper¬ 
vised  classification  technique  as  shown  in  Figure  1  and  Figure  2, 
respectively. 

To  obtain  a  better  result,  the  intersection  of  these  two  classification 
maps  (Figure  3)  was  used  as  the  final  decision  rule  for  defining  the  granite 
regions.  The  advantages  of  Figure  3  over  Figure  1  and  Figure  2  are  several: 

(1)  It  avoids  random  errors  in  the  decision  map  of  the  supervised 
classification  methods,  particularly  those  located  far  away 
from  cores  of  the  training  sets; 

(2)  It  has  information  for  labeling  segmented  regions  from  the 
unsupervised  classification  method;  and 

(3)  the  intersection  of  two  decision  maps  in  fact,  strengthens 
the  probability  of  correct  classification. 

In  this  report,  the  "intersection"  of  two  classification  methods  was 
done  by  use  of  an  overlay  analysis  of  two  decision  maps. 

Theoretically,  a  new  algorithm  should  be  developed  to  perform  the  task 
of  merging  two  image  processing  methods  centered  around: 

(1)  Class! fying  the  results  of  the  region  algorithm  using  the 
training  sets  information,  and 

(2)  Extracting  distinctive  regions  using  multiple  ratio  bands, 

and  classifying  them  according  to  the  training  sets  information. 

It  is  believed  that  certain  techniques  of  artificial  intelligence  are 


useful  to  this  effort. 
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2.  Fine-tuning  of  the  Region  Algorithm. 

a.  Dynamic  edge  patterns  as  indication  of  rock  types 

For  this  technical  report,  the  region  algorithm  was  utilized 
to  segment  terrain/rock  types  in  multiple  stages,  and  to  reveal 
the  evolutional  patterns  of  the  edges,  such  as  lineaments  or  con¬ 
tact  zones  between  two  lithological  types,  at  the  study  area. 

An  improvement  of  the  region  algorithm  can  be  made  by  adding 
a  subroutine  aimed  at  testing  whether  edges  within  a  larger  region 
are  noise  or  real  boundaries.  If  they  are  determined  as  noise,  they 
can  be  removed,  and  vice  versa.  This  capability  will  provide  the 
anlayst  with  a  sounder  basis  for  the  regionalization  of  terrain  and 
lithological  types. 

b .  Display  Texture  and  Tone  information  of  each  Segmented  Region 
To  provide  more  information  for  the  analyst  to  make  decisions 

regarding  identification  of  rock  types,  another  subroutine  can  be 
added  to  output  the  texture,  tone  and  size  information  of  the  seg¬ 
mented  regions  under  investigation.  Using  such  information,  the 
analyst  may  be  able  to  identify  the  terrain  and  rock  types  by  means 
of  a  comparative  analysis  provided  that  certain  texture-tone  charac¬ 
teristics  of  training  areas  are  known. 

In  practice,  the  analyst  would  use  both  edge  pattern  and  quanti¬ 
tative  texture-tone  information  for  rock  types  discrimination  and 
identification. 

c.  Generalization  by  a  Color  Prediction  Model 

Conventionally  color  monitors  are  used  as  a  device  for  generating 
color  composite  from  multi-channel  data  by  means  of  a  coding  process 


using  three  color  primaries,  red-green-blue  or  yellow-cyan-magenta. 
With  known  quantitative  texture  and  tone  information,  specific 
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colors  can  be  generated  for  a  given  rock  type  by  means  of  the  color 
theory.  For  instance,  if  there  are  two  parameters  (1  tone  and  1  tex¬ 
ture)  for  each  region,  a  specific  color  for  that  region  can  be  made 

by  assigning  the  tone  information  to  Red  domain  and  the  texture  infor¬ 
mation  to  Blue  versus  Green  domain.  If  three  parameters  are  available 
1  tone  and  2  textures,  the  specific  color  code  for  that  region  can  be 

generated  by  assigning  tone  to  Red,  texture  1  to  Green,  and  texture  2 

to  Blue. 

For  the  two-parameter  system,  different  color/tensity  codes  can 
be  generated  if  one  allows  the  tone  parameter  to  control  intensity 
levels  and  texture  to  control  colors.  In  an  8-bit  color  monitor 
system,  for  instance,  one  can  allow  the  tone  parameter  to  display  16 
intensity  levels,  and  the  texture  variable  to  give  16  different  colors 
The  combination  of  such  intensity  and  color  codes  should  allow  the 
researcher  to  predict  and  identify  certain  rock  types  using  given 
texture  and  tone  information  from  either  ground  truth  or  laboratory 
analysis,  or  a  combination  of  both. 
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APPENDIX  1 

Texture  and  Tone  Variables  of  the  Training  Sets  with 
the  Original  LANDSAT  MSS  Data 
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APPENDIX  2 


Texture  and  Tone  Variables  of  the  Training  Sets  with 


Ratio  Bands  of  LANDSAT  Data 
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